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We predict a resonant microwave absorption on collective vortex modes in a superclean d-wave 
superconductor at low temperatures. Energies of the collective modes are multiples of the distance 
between the exact quantum levels of bound states in the vortex core at lower temperatures and 
involve delocalized states for higher temperatures. The characteristic resonant frequency is larger 
than the cyclotron frequency u c but lower than the Caroli-deGennes-Matricon minigap A 2 /Ef; it 
has a y/~H dependence on the magnetic field and decreases down to oj c with increasing temperature. 
We calculate the vortex mass as a response to a slow acceleration. This mass is equal, b y the ord er 
of magnitude, to the mass of electrons inside the vortex "core" with dimensions £ ky £yj H&/ H ; it 
increases with temperature. We discuss the universal flux-flow regime predicted inu and show that 
it exists in a broader temperature range than it has been originally found. 
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I. INTRODUCTION 

The vortex dynamics in clean superconductors is deter- 
mined by localized states in vortex cores. Motion of a vor- 
tex excites transitions between the states; a competition 
between the relaxation rate 1 /t and the interlevel spacing 
controls the proportion between_dissipative and reactive 
forces experienced by the vortexH'B. In o?-wave supercon- 
ductors, the vortex dynamics is expected to be more intri- 
cate due to a peculiar structure of the vortex-core states. 
The presence of gap nodes introduces the most impor- 
tant difference in the structure of core states compared 
to an s-wave superconductor. As was shown in RefsHru, 
instead of, a well-defined quasiclassical Caroli-deGennes- 
Matricontj interlevel spacing u>q, the true quantum states 
in d-wave superconductors have a much smaller separa- 
tion between quantum levels, Eq, which depends on the 
magnetic field. As a result, there appears another param- 
eter Eqt which influences the vortex dynamics. Accord- 
ing to Ref.cl, a new universal regime can be reached in 
superclean superconductors with longitudinal and trans- 
verse components of the conductivity tensor independent 
of the relaxation time. It is realized when the relaxation 
rate is smaller than the average distance between the qua- 
siclassical energy levels (ujq) ~ T c 2 / Ep but larger than the 
separation between the true quantum-mechanical states 
E . 

In the present paper, we study both steady and oscilla- 
tory motion of vortices at low temperatures T <C T c and 
low magnetic fields H <C H C 2 in more detail using the mi- 
croscopic theory outlined in Refll. For a steady flux flow, 
we find that the universal regime is realized in a consider- 
ably broader region of temperatures and magnetic fields 
than what was originally predicted in ReflH. This con- 
clusion agrees with the results of Ref.El. The relevant 
parameter which governs the vortex dynamics is shown 



to be the energy of "collective modes" E co \. We calculate 
E co \ and show that it coincides with the true quantum 
mechanical interlevel spacing Eq if the energy of excita- 
tions is e <C T c y/ H I H C 2, and involves delocalized states 
for excitations with higher energies e 3> T C -\J H j 'H C 2- For 
an oscillatory motion of vortices which can be excited 
by a microwave irradiation, we predict a finite resonant 
absorption even for infinitely long relaxation time. The 
series of resonances occurs at multiples of E co \ with the 
resonant frequencies considerably lower than the average 
quasiclassical interlevel distance T^/Ef but higher than 
the cyclotron frequency uj c . They have a \H magnetic 
field dependence at very low temperatures. For higher 
temperatures, the series of resonances compresses down 
to an absorption edge which approaches uj c for T ~ T c . 

II. NONEQUILIBRIUM EXCITATIONS 

A. Spectrum 

We consider superconductors which have the coher- 
ence length £ much longer than the inverse Fermi mo- 
mentum, ppl; ^ 1. For these superconductors, a quasi- 
classical approximation is appropriate. A quasiclassical 
particle moves along a straight line parallel to its momen- 
tum which is conserved with the accuracy (pf£) _1 - The 
angular momentum is thus another approximately con- 
served quantity even if the order parameter does not have 
a cylindrical symmetry; it is now a continuous variable 
not directly related with labelling the quantum states in 
the vortex core. A quasiparticle passing at some distance 
near a vortex can become classically localized in the vor- 
tex core region. It will have an energy characterized by 
the momentum along the vortex axis p z , the momentum 
direction in the plane perpendicular to the vortex axis, 



1 



and by the impact parameter b = —/J,/p± coupled to the 
angular momentum p. 

For a d-wave superconductor without admixture of 
other components, the order parameter is 



A sin(2a)e 1 ' 



(1) 



where a is the angle between p± and the x axis which 
is taken along one of the gap nodes. The modulus of 
the order parameter Ao = Ao (/?,</>); at large distances, 
Ao = Aoo. We take the z axis along the vortex circula- 
tion, z = sign (e)H/H, and denote by vj_ the quasiparti- 
cle velocity in the (x, y) plane. If p and <j> are the distance 
and the azimuthal angle in the cylindrical frame, the im- 
pact parameter & of a quasiparticle moving through the 
point (p, cf>) and the distance s along the trajectory are 



b — p sin(0 — a) ; s — p cos(0 — a) . 



(2) 



The quasiparticle energy can be found using the quasi- 
classical Green function technique in the same wgy as it 
has been done for s-wave superconductors in Refill. Here 
we summarize the most important properties of the en- 
ergy spectrum of a quasiclassical particle in the vortex 
core. It is discussed in more detail in Appendix. The 
quasiclassical spectrum has several branches belonging 
to various values of the radial quantum number n. The 
anomalous branch with n — crosses zero of energy as 
a function of the impact parameter. For small impact 
parameters b <C £/| sin(2a)|, 



e (a,b) = 



2A^sin 2 (2a)L6 p±_buj c 



(3) 



where L = In [l/\ sin(2a)|] for b <C £ and L = 
ln[£v±/bv F \sm(2a)\] for f < b < £/| sm(2a)|; the sec- 
ond term is the magnetic energy where u) c — \e\H/m±c 
is the cyclotron frequency with mi = pj_/v±. Modulus 
of charge appears due to the choice of the z axis. Note, 
that the magnetic energy has the energy level spacingQ 
equal to the Larmor frequency lu c /2 rather than just to 
lu c as it would be the case for the Landau spectrum in the 
normal state. This is because the wave function is local- 
ized at distances of the general order of £ rather than at 
the magnetic radius pu ~ p±c/eH which would be the 
case in the normal state. Near the gap nodes, i.e., when 
sin(2a) is small, the low-lying states with energies much 
below the gap at infinity e <C Ao| sin 2a |, correspond to 
b <C £/| sin 2a | . The derivative of the energy with respect 
to the angular momentum, u>o(a) = pJ_ 1 [deo(a, b)/db], is 
the analogue of the Caroli-deGennes-Matricon minigapcl 
in an s-wave superconductor. 

For large impact parameters, b 3> £/|sin2a|, the en- 
ergy is 



A M |sin(2a)|-^ 



pib(jj c 



(4) 



eo(&) = sign (6) .^00,0^^1 - —r^ , ^ 

The particles are localized on lines passing through 
the vortex perpendicular to the momentum direction. 



T he localizati on length along the trajectory is s ~ 
VsVI S n (2«)| <C b. A similar spectrum was considered 
in Reffl. 

Eqs. (||) and (^) can be combined into a single inter- 
polation formula 



Ag sin 2 (2a)b 
|6|Aoo|sin(2a)| + /3v A 



p±buj c 



(5) 



where j3 is a constant. The fit to the large-6 region of 
energy gives j3 — 1/4; the best fit for low-6 region is 
reached with j3 — 1/(2L). We define an angle ao at 
which the particle energy e = eo(6) — p±bu> c /2 is equal to 
the gap e = A^l sin2ao|; the angle ao is close to one of 
the nodes if e <C A m . It is a quasiparticle moving at an 
angle a > ao which is classically localized. 

Other branches with nonzero radial quantum numbers 
n are separated from eo(b) by energies of the order of 
Aoo | sin(2a)|. The number of branches with energies be- 
tween e and A^j sin(2a)| is 



j sin(2a) 



[A*, sin 2 (2a) — I 2 ] 



(6) 



Eqs. (^) and (|5|) are quantitatively correct as long as 
b <C Pmax where p max ~ 6 \] H C 2 / H is the distance be- 
tween vortices. Such distances are, in principle, accessi- 
ble for a particle which move s near th e direction of a gap 
node within the angle a ~ y/H/H C 2- We shall see, how- 
ever, that the condition b <C /3 m ax is fulfilled for energies 
e^A^y/H/Hcz. 

A classically localized particle is not necessarily local- 
ized in a strict quantum-mechanical sense. This was also 
found in numerical calculations of Ref.EI. The true quan- 
tum levels can be obtained frap 1| the semiclassical Bohr- 
Sommerfeld quantization rulcoEJiij 



p(a) da — 2n ( m + 



(7) 



m is an integer, | appears because the single-particle 
wave function changes its sign after encircling a single- 
quantum vortex. The angular momentum p(a) is ex- 
pressed through the quasiparticle energy according to Eq. 
(||). Consider first an energy e <C A^^J Hj H c i. The in- 
tegral in Eq. (0) converges according to Eq. (||) and 
is determined by angles a ~ y/H/H C 2. For this region 
of angles, the spectrum is actually given by the exact 
equation (||): eo(p) = — u>o{a)pL where 



cuo(a) = 8iloa 2 -I- ui c /2 



with 



O = (A^Mjjj.) H(A 00 /e) x /H/H c2 ] 



(8) 



(9) 



The characteristic impact parameters are of the order of 
b ~ e/pxLJc, i.e., b < p ma x. We obtain E m = E (m+l/2) 
where 
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E = \/n 



(10) 



B. Balance of forces on a moving vortex 



A qualitative expression for the spectrum with the v H 
dependence on the magnetic field similar to Eq. jl0| ) 
was found in Ref.Ea neglecting the magnetic energy. The 
argumentation was that the divergence ofu>Q 1 (a) with- 
out the magnetic energy should be cut off at such an- 
gles a that the characteristic distance of a qusiparticle 
from the vortex axis is of the order of the intervortex 
distance p max ; here the corrections to the quasiparticle 
energy induced by neighboring vortices is just of the or- 
der of u>o(a). We see, however, that Eq. ( |i~0| ) is obtained 
exactly if one takes into account the magnetic quantiza- 
tion. We find simultaneously that the particles with en- 
ergies e <C Ay/ H/H C 2 are acually localized in the vortex 
core since b <C p m ax! thus the corrections from neigh- 
boring vortices arc not important. For larger energies 
e ~ Aqo y/H/Hcz, the localization radius becomes of the 
order of p max . These particles can not be considered as 
localized any more because of the presence of other vor- 
tices. For higher energies e 3> A ooV / i//iJ C 2, the part 
with extended trajectories a < ao in Eq. (Q) dominates, 
and thus quasiparticles are no longer localized in the core 
even for an isolated vortex. 

The average quasiclassical energy-level spacing (uj a ) ~ 
flo ~ A^/Ep determines the parameter Q q t which con- 
trols the vortex dynamics in clean s-wave superconduc- 
tors with t T~ x . The value CIqt ~ 1 separates the 
so called moderately clean limit, fioT <C 1, with a highly 
dissipative vortex dynamics from the superclean limit, 
f2o T 3> lj where dissipation is small. 

In d-wave superconductors, ^situation is more intrigu- 
ing. It was predicted in Ref.El that a large dissipation 
persists for much longer I. In the present paper we show 
that, for a steady vortex motion, the relevant parame- 
ter which marks the transition between dissipative and 
nondissipative dynamics is E co \t where E co \ is the tem- 
perature dependent characteristic energy of "collective 
modes" induced by moving vortex. For low tempera- 
tures, it coincides with the energy of "bound states" de- 
termined by Eqs. (@), @: E col = E < O . The 
energy E co \ decreases down to w c when the temperature 
approaches T c . The parameter uj c t is much smaller than 
Hqt since lo c ~ (H/ H C 2)£Iq. We have thus a hierarchy of 
energies uj c <C E co \ <C SIq. To get the parameter which 
controls the vortex dynamics, these energies should be 
compared either with the relaxation rate 1/r for a steady 
motion of vortices or with max (1/r, to) where u> is the 
characteristic frequency of vortex oscillations. Large val- 
ues of £Iqt require a very high purity of samples: the 
mean free path should be £ 3> (T c /Ep)£. Nevertheless, 
there are experimcntal.-c.vidences that such a regime can 
be realized in practice£3 



The force on a vortex from the environment where it 
moves is the momentum transferred from the excitations. 
The exact microscopic expression for the force has been 
derived in Ref.lj. It can equivalently be presented as a 
force produced by quasiclassical particles with a distri- 
bution /, characterized by canonically conjugated coor- 
dinate a and angular momentum The contribution 
from classically localized states is 



p(loc) 



E 



dp z dfj, da dpn 
~2tt 2tt ~dt 



(11) 



Using dp/dt = — Ve„ = [z x p](de n /d/i) we get 



p(ioc) _ 



dp z dfi da 




z x p±\—df n 



■kN[z x u] 



-irNu 



dfo de n 
de d/i 



dfo de n 
de a/j, , 



d/j, 



(12) 



We take the nonequilibrium distribution function in the 
form 



5/ = -d7 ((upj 



)7H + ([u x p_l]z)7o) 



(13) 



where fo is the equilibrium Fermi distribution, the fac- 
tors 70 and "/h describe the longitudinal and transverse 
responses to the vortex velocity u, and N is the electron 
density; (• • -) F is the average over the whole Fermi sur- 
face which includes averaging over a. We assume here a 
Fermi surface with not less than the tetragonal symmetry 
in the plane perpendicular to the vortex axis. For sim- 
plicity, we consider only the particle-like Fermi surface, 
the generalization for a surface with both particle-like 
and hole-like parts is given later. n 
The force from classically delocalized states isa 



p(del) _ 



dfo 



= ttN[z x u] ( / ^h-^t de 

Vdel de 



-nNu 



/ ^fo-j-de 
Jdd de 



(14) 



The force F onv should be balanced by the Lorentz force 
Fl = sign(e)0o[j x z]/c where 4>o = irc/\e\ is the flux 
quantum. Since the average electric field is E = [B x 
u]/c, the force balance F^ + F cnv = gives a linear 
relation between the transport current and the electric 
field. The proportionality coefficients are the Ohmic and 
Hall conductivities. The Ohmic conductivity is 




dfo de n 
r Yo- r -z—dfi 
de Ofi i 



dfo, 
70— de 



de 



(15) 
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For the Hall conductivity an one has the same expression 
with jh instead of 70 and \e\ replaced with e. 

In Galilean invariant systems j = N s ev s , and the force 
balance is sometimes presented in terms of the Magnus 
force from the superfluid component and the friction plus 
transverse forces from the normal component 

ttN s [(v s - u) x z] - Du- D'[i x u] = 0. 

In such representation, the Magnus force includes the 
Lorentz force and a part of the transverse force F onVj j_. 
The constants D and D' are expressed through the con- 
ductivities as 

c 2 D . ( .c\nN s ~D') 

4> b 4> b 



III. DISTRIBUTION FUNCTION 

A. Kinetic equation 

The kinetic equation for the distribution function 
f(t,a, (m) for a system of fermions characterized by-canon- 
ically conjugated variables fi and a has the formed 



dj_ 

at 



Of dt n 

da d\i 



de n df 
da dfi 



f-fo 



(16) 



A particle is classically localized if it has an energy e < 
Aqo and moves at an angle a > a$ counted from one 
of the nodes. For localized quasiparticles, the derivative 
de n jd[i = —oj n (a) is the interlevel spacing. The collision 
integral is written in a r-approxima,tion. This equation 
can also be derived microscopically.0 

If the vortex moves with a velocity u with respect 
to the heat bath, the Doppler shift of the energy e — > 
e(fi, a) — pu produces the "driving force" (de n /da) = 
[u x p_l]z acting on quasiparticles. The third term in the 
l.h.s. of kinetic equation ( |l6| ) contains this force multi- 
plied by dfo/dfi = —uj n (a)(dfo/de). We look for a so- 
lution in the form / — /q = Sf where Sf — 5f(e,a) is 
independent of b. Eq. (|l^) gives 



dSf 

3a-- ([UX 



pjz)^ = U(a)Sf 



where 



(17) 



(18) 



Here we assume that the vortex velocity has a form u = 
u u e~ luJt where u <C A. 

Introducing the longitudinal and transverse responses 
to the vortex velocity according to Eq. (|l3|) and tak- 
ing into account that, for the tetragonal symmetry, the 
responses do not depend on the direction of the vortex 
motion with respect to the crystal lattice, one finds two 



coupled first-order differential equations for 70 (a) and 
1h(ol): 



die 



-JH-U (a)jo + 1 = 0, 



da 



+ 7o - U(a)iH = 0. 



(19) 



A quasiparticle is delocalized either for energies e > 
Aqo or for angles a < a$. For a homogeneous magnetic 
field the distribution function of delocalized electrons was 
shownp to satisfy Eq. ( |l7| ) where u) n is replaced with gu> c . 
Here g is the number of states for a particle with given 
e, a, and b at large distances from the vortex, i.e., the 
quasiclassical Green function g: 



g = l/\Jl 2 - Ag sin 2 2a. 



The distribution function thus has the form of Eq. (|1 
with 70 and jh satisfying Eqs. ( |l9| ) where now 



U(a) = ( -iw + - ) {u c g)~ 



(20) 



Since 70 and jh obey first-order differential equations 
they are continuous functions at a — uq. For e > A^, 
the potential is given by Eq. (^0|) in the whole region of 
angles. 



B. Static response 

Equations ( pj| ) can be easily solved. We have 7h(«) = 
Re IF; r fo( a ) = ImVF where 



W = e [ iQ+F ( Q )i [ C 



'f 

Jo 



[ M ' +F ( Q ')] da' 



with 



F(a) = / U{a')da. 
Jo 

In the moderately clean limit, VIqt <C 1, the potential 
U{a) is always large, and we_pbtain the local solution as 
in an s-wave superconductora 



70(a) = uj {a)T, 7#(a) = [w (a)r] 



(21) 



In the superclean limit Qqt <^ 1, the potential U(a) 
is small almost everywhere except for vicinities of the 
gap nodes where it becomes large. Consider first ener- 
gies e <C Aqo yj H / H C 2 and find the distribution function 
for the anomalous branch n = in the region of angles a 
not specifically close to the gap nodes. It is this branch 
which is only excited at temperatures T <C T C \J H/H C 2- 
The overall behavior of the distribution function is to the 
highest extent determined by what happens in a close 
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vicinity of the gap nodes. It is this region which is re- 
sponsible for the whole build-up of nonequilibrium dis- 
tribution of excitations. 

For angles larger than a ~ \J H/ H c2 from the nodes 
one can neglect the potentia l U. As a result, one has for 
Sa < a < 7r/2— 5a where y/H/H c2 <C 5a <C 1, 



70 (o0 = ^4 cos a + B sin a; 
7fj(a) = 1 — A sin a + Bcosa. 



(22) 



The constants A and -B can be found by matching with 
the solution in the vicinity of the node, a <C 1, where 
1o,h{o) oc e f '°'. This provides the boundary condition 
1o,h{+5o) — e 2X jo,H{—5a) across the node at a = 0. 
Here 2A = F(+5a) — F(—5a). For such energies, the 
region of angles a < ao with delocalized trajectories is 
not important. The integral for A converges and is de- 
termined by angles a ~ yj H/H c2 . This range of angles 
sets the width of the transition region near a gap node 
where the distribution function jumps from its value at 
a = —5a to its value at a — +5a. We obtain using Eq. 

(i 



A = r _1 



da 



4Enr 



(23) 



The solution for jo,h should be periodically contin- 
ued to the rest of angles with the period tt/2 since the 
response function has the same tetragonal symmetry as 
the underlying system: j(tt/2 — 5a) = j(—5a). Together 
with the above boundary condition, this gives 



A = 



e x sinh A 



,-A 



2 sinh 2 A + 1 



B = 



sinh A 



2 sinh A + 1 



(24) 



We have after averaging over the azimuthal angle a 



(in) 



4 tanh 2 A 
7r 1 + tanh 2 A 



. (7o) 



4 tanh A 
7r 1 + tanh 2 A ' 



(25) 



For energies e 3> A^y/ H/H c2 , the region of angles 
a < ao with extended trajectories dominates; now it is 
ao which determines the width of the transition region 
where 70 and jump. Compact expressions can be 
obtained if we assume that U(a) is independent of a for 
a < ao. For simplicity, we replace guo c T with uj c t. This 
does not change the results qualitatively. We have for 
a < ao 



70 = di + e a/u)cT (Ci sin a + C 2 cos a) , 
1h = d 2 + e a ' UcT (Ci cos a-C 2 sin a) 



(26) 



where 



di = 



LO c T 



1 



rf 2 = 



1 



We can use Eq. (|22|) for angles a > ao + 5a. Since the 
contribution to A form the region of angles a ~ 5a is 



much smaller than that from the angles a ~ ao, we can 
extend Eq. (53) to the angles a = ao and a = tt/2 — 
ao. First, we match equations ( |22] ) and ( p6| ) at a = an. 
Another condition is obtained by matching Eqs. (|2q ) 
taken at a — —ao with Eq. ( p2| ) taken at a = tt/2 — ao- 
We obtain four equations for four constants A, B, C\, 
and C2 where now A = ao /lu c t — e/(2A 00 cj c r). 

If A is not considerably smaller than unity, we can ne- 
glect ao in the boundary conditions, and obtain 



.4 : 
B - 



(1 
(1 



sinh A 
2 sinh 2 A + 1 



-X 



-A 



sinh A 



2 sinh 2 A + 1 



di 



di 



sinh A 



2 sinh 2 A + 1 

e A sinh A 
2 sinh 2 A + 1 ' 



For values of A ~ 1 and larger, one automatically has 
lu c t < 1. As a result, both d\ and d 2 are small, and we 
recover Eqs. (^4j) and ( p5| ) with the new definition for A. 
The coefficients G\ and C 2 are in this case 



Ci = 



cosh A 



2 sinh 2 A + 1 



C 2 = 



sinh A 



2 sinh 2 A + 1 



We s ee that the distribution function for energies 
Aooy/ H / H c2 <C e Aoo is not small even for delocalized 
trajectories a < ao- 

The contribution from the region of angles a ~ 1 to 
Eq. ( |l5|) and to the corresponding equation for <th de- 
creases together with A as the parameter oj c t increases. 
For A -C 1, the contribution from angles a ~ ao where 
the gap Aoo| sin(2a)| is of the order of temperature be- 
comes important. In this case, we obtain from Eqs. ([22]) 
and @ 

A = A(l - d 2 - di) + a (l -d 2 + di); 
B = A(l — d 2 + di) - a (l - d 2 - di); 
C 1 = l-d 2 + Xd 1 ; C 2 = A(l - dj,) - d x . 

As a result, 
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(70) = -A; (7^) = 1. 



(27) 



We can thus use Eq. (|25|) within the whole range of A 
and combine the two results for A in different regions of 
energy into a single interpolation expression 



A = 



(28) 



where E co \ has the meaning of a characteristic energy of 
"collective modes" 



1 



2kl 



which we discuss later. 

For energies e > A M we have from Eqs. (|2 
C 2 = and 

7o = di;7H = d 2 . 



(29) 

Ct = 
(30) 
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C. Frequency-dependent response 

Assume that the external frequency lu <C f2o- I n this 
limit, the solution of the kinetic equation ([r]) has the 
form of Eqs. ( p5|p8| ) where 1/t in Eq. ( p8| ) is replaced 
with 1/t — ilu. Consider first the case u> 3> lo c . The 
factors d\ ~ lo c /oj and d 2 ~ w 2 /w 2 are small. In the 
limit wt — > oo, the longitudinal response is 



(7o) = 



. 4 tan A' 
'tt tan 2 A' - 1 



where 



A' = 



4£ C oi 



(31) 



(32) 



The response has poles at A' = ir/4 + irk/ 2 or at fre- 
quencies lo = LUk where lj^ — E co \(l + 2k) and E co \ 
is given by Eq. (^9|). Harmonics with k 7^ appear 
due to the absence of the axial symmetry. N ote that 
E co \ ^> uj c . For low temperatures T/T c -C yjHjH c2l 
the eigen frequencies are independent of temperature: 
ujk = Eo(l + 2k). These poles are the collective modes 
of electrons involved into the vortex motion. For low 
energies e <C A^^J H/ H c2 , these modes coincide with 
multiples of the distance between the true quantum me- 
chanical energy levels determined by Eq. (|7]) . For higher 
energies, delocalized quasiparticles are involved into the 
vortex motion. 



In the superclean limit flo T 3> 1 the factors 70 and 
7# for low temperatures T <§C T c are determined by Eq. 
(p5|). The general expression for the conductivities are 
rather complicated. We consider two limits. The uni- 
versal regime is reached when A 3> 1, i.e., when cither 
^H/T172 < 1/ (fior) for T/T c < ^H/H^ or when 
lj c t < T/T c for T/T c > ^/H/H c2 . The region of the 
universal regime is thus much larger than it was pre- 
dicted in Ref.II This was pointed out by Makhliro. The 
condition A > 1 automatically implies that uj c t <C 1, 
therefore, both d\ and d 2 are small. One has from Eq. 
( p5| ) ju = 1 -p.2/71", 70 = 2/ ir which results in universal 
conductivities!!! 



*0 = Jf(*e 



TT 

N h ) ( 1 - - 



In the limit A « 1 we have from Eq. fl27| ) 

(AT e - iV/j) ec 



ec(N e + N h ) A 

o"o = -p. ; &h 

B T 



B 



(34) 



(35) 



where A = {N e K e + N h A h ) (N e + N h ) 1 and 

cosh- 2 ^ de 
4Tln2 



A e ,ft = 





/ 1 



(E C ol) Fe, I'h ' " \ 2T 1 'IT 



(36) 



IV. CONDUCTIVITIES 



A. Steady motion 



Since 70 and 7// do not depend on fi, the integration 
over d/i in Eq. (15) can be reduced to the integration 
over de. In the sums over n, only the term with n = 
remains because all w n with n / are odd functions of 
\i. If the Fermi surface has electron- like and hole- like 
pockets, the conductivities are (compare with Ref.cl) 



00 



|e|c 
B 
ec 



F,e 



dc 

dfo 
de 



de- N h 



(1h) 



F,h 



de 

dfo 
de 



de 



In the moderately clean case, the factors 70 and 
are given by Eq. (Ell) . The conductivities are 



B E F \T J \2T 



E F 



HTc/T) 



(33) 



where N ~ A e + iV^ 
s-wave superconducto: 



similar to the results for an 



Here J7o is determined by Eq. (jRj) with the logarithmic 
factor \-a[{T c /T)y/H/H c2 ]. The average is taken over the 
particle-like and holejhke parts of the Fermi surface. Eq. 
(|35| ) also agrees withO. This regime is reached when ei- 
ther yjH/H c2 > 1/ (fi T) for T/T c < ^H/H c2 or when 
w c r » T/T c for T/T c » ^/H/H^. 

As the temperature approaches T c , the contribution 
from fully delocalized states with e > A^ becomes more 
and more important. The delocalized states give the 
normal-state Ohmic and Hall conductivities in the limit 
T — > T c . In the limit lu c t <C 1, the flux-flow parts of 
ao and an start from the universal values and then first 
decrease as A^/Tc ~ y/l — T/T c as long as f2 r remains 
large. The ratio A^/Tc measures the relative contribu- 
tion from localized states. With T further approaching 
T c , the moderately clean regime is reached, and ao be- 
comes proportional to (A 00 /T c )0o oc (1 — T/T c ) 3 ^ 2 while 
a H is proportional to (A^/Tyfi 2 , cx (1 ~^/T c f /2 being 
small: an /cro ~ ^o T (compare with Ref .□) . In the limit 
lo c t 3> 1, the universal regime is not realized, and the 
conductivity ao starts from Eq. (|35|) and saturates at 
ao ~ {Nec/ B){oj c t)~ 1 while Off remains constant. 



B. Dispersion of conductivity 



G 



1. Low frequency: vortex mass 



2. High frequency: resonances. 



If both the frequency and relaxation rate are low such 
that A, A' <C 1 the transverse response is ("fn) = 1- The 
transverse force is 



F env ,± = 7T (iV e - N h ) [z x u 
The longitudinal response, 
4 



(37) 



<7o) = 



irE, 



col 



gives the force 

Fonv.n = -Du 
at 

where D is the friction coefficient and M plays the role 
of the vortex "mass" . Indeed, the force balance takes the 
from of the Newton's law 



- Du = M 



r du 

w 



(38) 




Here the vortex acceleration is a result of action of the 
Lorentz, transverse, and friction forces. For a steady flux 
flow without a friction, we would have from Eq. (|3^ ) 
j s = [N e — Nh) eu. In a system with the Galilean invari- 
ance, where N = N e and Nh = 0, the sum Fl + F onv ,j_ = 
iJV[(v s — u) x z] is the Magnus force. The friction coef- 
ficient and the mass per unit length are 

D = 7T (N e + N h ) A/t; M = tt (N e + N h ) A. 

The effective mass of a d-wave vortex 



M ~ Nmi 2 



is much larger than the mass of a conventional vortexEl 
M ~ Nm^ 2 obtained in the same-way. For temperatures 
T/T c < ^H/H c2 , it is the masai3 of electrons inside the 
"vortex core" with the dimensions £ by The 
mass increases with temperature. Note that this mass 
is very large compared with what is usually obtained 
by calculating the corrections to the kinetic energy of 
a moving vortex proportional to the second power of its 
velocity including the electromagnetic energy. For exam- 
ple, the recent result of Ref.EZl is roughly by the factor 
of (T c /Ef) 4 smaller! We have to stress that this mass 
appears as a response to a slow acceleration. With an 
increase in the characteristic frequency of the vortex mo- 
tion, the response transforms into a highly dissipative 
resonant behavior where a mass is not an appropriate 
quantity. 



Due to the poles in (70), the dissipative component 
of the response has a real part even for r — > 00. For 
T <C A 00 y / oj c /H,o, the absorption is concentrated near 
Uk{Pz) = (2fc + 1)Eq{p z ) with the line width Sujk/^k ~ 
(T/AooVftoM < 1: 



Re ao 



Nec^ /4E Q 



B 



Averaging over the Fermi surface, i.e., integration over 
dp z , gives rise to the van Hove singularities in conduc- 
tivity at frequencies where dEo/dp z = 0, i.e., where 
lo = {2k + l)E { P z =0). 

For T 3> Aqo \J ^c/^Oi we have 



Recro 



Nec / A. 



B 



Tuj 



— cosh 

fe>0 



7tA p 



(1 + 2*) 



The conductivity is exponentially small for low frequency 
uj -C Eg = ttAooU! c /AT. It starts to rise at the absorption 



edge w ~ E g . The average in the r.h.s. of Eq. (39) is 



equal to 2 /it for large frequencies w » £ s , when transi- 
tions between many resonance modes are excited. This 
is equivalent to the limit A > 1 for a steady flow. 

The absorption occurs at the collective modes of elec- 
trons involved into the vortex motion. For temperatures 
T T C ^H/H C 2, these modes are transitions between 
the true quantum mechanical energy levels determined 
by Eq. (M). For higher temperatures, delocalized quasi- 
particles are involved into the vortex motion, and the 
resonances take place at collective modes with energies 
of the order of E„. 



V. SUMMARY 

For a steady vortex motion, one can identify three 
regimes in the order of increasing r: (i) The moder- 
ately clean regime with £Iqt <C 1. In this case, the 
main response of vortices to the external current is dis- 
sipative with only a small Hall angle proportional to 
f2o T - The conductivities behave similarly to the s-wave 
case, see Eq. ([33j). Next comes the superclean limit 
floT 3> 1. However, in d-wave superconductors, it fur- 
ther separates into two sub-regimes. The relevant pa- 
rameter is the energy of collective modes excited by the 
moving vortex. The vortex dynamics in a cZ-wave su- 
perconductor depends crucially on the behavior of ex- 
citations near the gap nodes. Due to the presence of 
gap nodes, only the excitations with very low energies 
e <C Aqo y/H/H C 2 are localized in a vortex core. The ex- 
citations with higher energies, coming into play for tem- 
peratures above T c y/H/H c2 , are actually the collective 



7 



modes where both classically localized and delocalized 
particles participate. The interplay between the energy 
of these collective modes E co \ and the relaxation rate or 
the frequency of vortex oscillations determines which of 
the two sub-limits is realized: (ii) Intermediate or "uni- 
versal" regime with E co \t -C 1 but Q,qt ^> 1. Here both 
dissipative and Hall components of conductivity tensor 
are universal. They are independent of the relaxation 
time and are only determined by the magnetic field and 
by the number of electrons and holes under the Fermi sur- 
face in the normal state, Eq. (|34j). (iii) Extremely clean 
limit E co \t 3> 1 where the dissipative part of the vortex 
response is small, and the transverse Hall component of 
conductivity dominates. 

Provided the superclean condition Qqt 3> 1 is satis- 
fied, the universal regime can be realized for tempera- 
tures T/T c <C y/ H/ H c2 under the condition E t <C 1, 
i.e., f2 T <C y/H c2 /H. This condition for the universal 
regime .was predicted in Ref.H However, as was noticed 
in Ref.S, the universal regime is not restricted to this 
temperature range. It can also be realized for higher 
temperatures if oj c t <C T/T c . For these temperatures, 
the dissipative conductivity vanishes only in the extreme 
clean limit lu c t 3> T/T c . 

The temperature and magnetic-field dependences of 
conductivities in the extremely clean limit are quite un- 
usual. Extremely clean limit can be reached by increasing 
the magnetic field above H/H c2 ~ (£!ot)~ 2 for T/T c < 
(flor)- 1 or above H/H c2 ~ (O t) _1 (T /T c ) for tempera- 
tures (f2o T )~ 1 < r/T c < 1. At the crossover from uni- 
versal to extremely clean regime, the Hall conductivity 
grows by the factor ir/(ir — 2) w 2.75 while the Ohmic 
conductivity decreases, Eq. (|35|). Fori temperatures 
T/T c < y/H/H c2 ^ its field dependence isi a oc #~ 3/2 ; 
it transforms into! a oc iJ~ 2 for T/T c > yfH/H c2 . The 
Hall angle approaches tt/2. 

At temperatures T ~ T c but not specifically close to 
T c , both Ohmic and Hall conductivities are of the or- 
der of Nec/B as long as lu c t -C 1. They, of course, do 
not have the universal values any more because all an- 
gles contribute similarly to the conductivities. However, 
the main conclusion is qualitatively the same: dissipa- 
tion only disappears when lj c t 3> 1. This is the major 
difference between s-wave and d-wave superconductors: 
in former, dissipation vanishes already for £Iqt 3> 1. For 
T — > T c , the superclean condition A^r/Ep ^> 1 fails, 
and one approaches the moderately clean limit where dis- 
sipation dominates. 

Having excited an oscillatory motion of vortices by a 
microwave irradiation at temperatures T/T c <C y H/H c2 
with the frequency u> ~ (A 2 ^/ Ep)-J H/ H c2 , one can ob- 
serve absorption resonances at the vortex-core states, if 
the condition lot ^ 1 is satisfied. The required frequency 
is higher than the cyclotron frequency and has a yff de- 
pendence on the magnetic field. For higher temperatures 
T/T c ^s> \/H/H c2 , one observes an absorption edge at a 
frequency of the order of (T c /T)uj c where the absorption 



increases sharply with the increasing frequency due to 
resonances at the vortex collective modes. 

Finally, we find that the response of a vortex to a slow 
acceleration is characterized by a "mass" term. The mass 
per unit length, by the order of magnitude, is equal to the 
mass of electrons inside the vortex core with dimensions 
£ by £y H c2 /H and increases with temperature, i.e., it is 
enormously larger than the electromagnetic contribution 
or any other correction to the kinetic energy of a vortex 
proportional to the second power of its velocity. 



ACKNOWLEDGMENTS 

I am grateful to Yu. G. Makhlin and G. E. Volovik 
for illuminating discussions. This work was supported 
the Swiss National Foundation cooperation grant 7SUP 
J048531 and by the INTAS grant 96-0610. The support 
by the Russian Foundation for Fundamental Research 
grant No. 96-02-16072 and by the program "Statistical 
Physics" of the Ministry of Science of Russia is also ac- 
knowledged. 



APPENDIX A: QUASICLASSICAL SPECTRUM 
IN D-WAVE SUPERCONDUCTORS 



The quasiclassical approximation, a quasiparticle pass- 
ing near a vortex has a definite trajectory which is actu- 
ally a straight line characterized by an impact parameter 
b. This argumentation holds for any symmetry of the or- 
der parameter and can be applied to a d -wave supercon- 
ductor, as well. We solve the quasiclassical Eilenberger 
equations for the Green-Junctions g, f , and using the 
same scheme as in Rcf.13. In the frame with the impact 
parameter b and the distance along the t rajectory s as the 
coordinates introduced in Section II A, the Eilenberger 
equations are 



df 

-iv±-j± -2e7 + 2A p .g = ; 
os 

iv F ?J--2ep + 2Alg = . 
os p 



(Al) 
(A2) 



and g 2 — f f< = 1. Here the functions g, /, and f> 
are either retarded or advanced Green's functions. The 
energy e = e+ (2e/c)v±A s . For an extreme type II super- 
conductor with the Ginzburg-Landau parameter k> 1, 
the magnetic field is nearly homogeneous, H = n£0o 
with = Hp/2. Since the vector potential component 
A s = —{b/p)A ( f > we have e = e — eHbv±/2c. 

The order parameter of a <i-wave superconductor has 
the form of Eq. (h]) We introduce 



V = ~Jo ~ds~ dS= 2 +a ~^( s )- ( A3 ) 
so that cos 7/ = b/p and sin?] = s/p. We put 
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f = fo exp(i(f) + irj); p = fl exp(-i0 - irj), 



and 



/o = -*(*)+ *C(«); /o t = ^) + iC(«) 



The functions and C satisfy the equations 
d( 

v±— — h 2e8 — 2i Ao P g sin n = 0; 
as 

— 2eC — 2iAo P g cos 77 = 0. 

as 



(A4) 
(A5) 

(A6) 
(A7) 



where g 2 = 1 — Q 2 — Q 2 . 

For energies close to the eigen energy, i.e., near the 
poles, the Green functions / and g are large. The func- 
tion g should be proportional to £ which has the required 
symmetry with respect to s. Therefore, we assume that 



C 2 > 1 



(A8) 



so that g — i^J C, 2 + 6 2 — 1 w i(. The sign is chosen to 
ensure a decay of g at large distances from the vortex. 
The function £, being the amplitude of the residue of g 
for the pole, decreases and should completely vanish at 
large distances.. It is reasonable to assume that, at large 
distances, £ — > and 9 2 — > 1. At large distances from 
the vortex axis, the full Green function 



f R = 



A r 



VIA P | 2 -e~ 2 



(A9) 



In combination with Eq. (A5), it gives the required sign 
of 8. Finally, the boundary condition for the Eilenberger 
equations are 



C — > 0; 6 — > ±sign [sin(2a)] for s — > ±oo. 
The solution is 

,R(A) = sign [sm{2a)]v_ L e- K 

X^oo — sm (2a)| cost? ± iS] e~ K ds 

9 R W = 2sign [sin(2a)] 

/ A |sin(2a)|cosr/] e - K ds' 



(A10) 



where 

A" 

As a result, 



J_ [e — Ao| sin(2a)| cos 77 ± i<5] e _if ds 2 / 
A |sin(2a)| shirks'. (All) 



iv ±e 



-K 



where 



C* [e - e ± i<5] 
A 

-OO 



(A12) 



cos rje 



-K 



ds 



eHbv ± /2c 



(A13) 



with 



C 



e~ K ds. 



(A14) 



Eq. ( |A13| ) was obtained ini-0 for an s-wave superconduc- 
tor. 

For small impact parameters, b <C £, one has | s |= p. 
Therefore, 



eo(&) = M C 



- l |sin(2a)| / 
Jo 



A n e-^ ^ 



\ e " ' hp±w c /2 



(A15) 



where oj c — \e\H/m±c is the cyclotron frequency with 
mj_ = Modulus of cha rge appears duetto the 

choice of the z axis. Eq. (A15) was obtained iifl. Due 
to the Kramer and Pesch effect^, the core size shrinks at 
low temperatures to £i ~ (T/T c )£. The spectrum takes 
the form of Eq. (g) with L = In [l/\ sin(2a)|]. 

The contribution from the magnetic field increases if 
the number of vortices increases together with the mag- 
netic field. However, the order-parameter phase gradient 
V% at the position of the vortex core under considera- 
tion remains equal to that produced by a single vortex 
^4>X = 1/p as long as p < 

pmax where p m ax is the in- 
tervortex distance since contributions from other vortices 
cancel each other. 

Near the gap nodes, where sin(2a) <C 1, Eq. ( |A13| ) 
for the anomalous branch can be simplified. Significant 
values of K are determined by p 3> £. Therefore, one can 
put A = AoqI 



K{s) 



2Aool sin(2a) 



2A r 



| sin(2a) 



s' ds' 



Vj_ 

The normalization constant is 
yp/b _pdp_ 



(p-b) 



C = 2ei 



Vp^p 



= be<K 1 { 1 ) (A16) 



where K\ is the Bessel function of an imaginary argument 
with 7 = 2AOO&I sin(2a)|/v^. The energy becomes 



eo(b) — A | sin 2a 



go(7) 
^1(7) 



■ bp±uj c /2. 



(A17) 



As a result, for £ -C b -C £/| sin 2a |, we get Eq. (||) with 
L = \n[^vj_/bvp\ sin(2a)|]. The low-lying states with en- 
ergies much below the gap at infinity, e -C Ao|sin2a| 
correspond to b <C £/|sin2a|. For large impact param- 
eters, b 3> ^/| sin 2a I, we can use the asymptotics of the 
Bessel f uncti ons for large 7 and arrive at Eq. (|4|). 

Eq. ( |A17| ) holds when Eq. (|Aj) is fulfilled which is 
equivalent to 



(e — A cos rf)e K ds <C v± 



(A18) 







for such s that K ~ 1. Eq. (A18) is satisfied either for 
7 C 1 or for 7 > 1. In the l atter case, t he exponent 
is localized at distances s ~ \/bj/\ sin(2a )| <C b. Using 
Eq. (H) we find that the integral in Eq. ( A18 ) is of the 
order of v±s/b ~ v± \/S,/b\ sin(2a)| <C v±. We see that 
for b 3> £/|sin(2a)|, a particle is localized at cos 7/ = 1. 
i.e., on a line passing through the vortex perpendicular 
to the momentum direction. 

The energy levels for n ^ can be found using the full 
quasiclassical scheme. It works either for b ^> £ or for 
n ^> 1. We look for a solution for Bogoliubov wave func- 
tions in the form e l I pdp . The quasiclassical momentum 




bv±\ 



< + vJ ~ A ° sin2(2a) (A19) 



where q p = p± yl — b 2 / p 2 . The integral along the closed 
quasiparticle trajectory is 



The separation between the gap and e increases with b 
until Xq ~ 1. The derivative becomes positive for larger 
impact parameters when po — > 6, i.e., when -C 1. The 
energy is 

e„ = Aool sin(2a)| - — . 

The correction to A is two times larger than that ob- 
tained for eo(b) from the more exact Eq. This is 
because the full quasiclassical approach does not work 
well near the turning point. The minimum of energy is 
reached at xq ~ 1 and is equal to 



1 



A 00 |sin(2a)| 1 - 



2ir 2 n 2 



The number of branches with the energy between e and 
Aqo is thus given by Eq. (^|). 



pdp 




2+0) - A^ur(2,.ir/ /) 



(A20) 



Consider positive and large t > ^ and positive ener- 
gies e — > Aoo| sin(2a)|. At large distances, 1 — Aq/ Aqq ~ 
£ 2 1 p 2 . For b > i ; , the term with 6 in Eq. ( |A20| ) domi- 
nates over the large-distance correction to Aoq. The zero 
of the square root in p determines the turning point 



bv± 



2(A 00 |sin(2a)| - e) ' 



For po b , the integral in Eq. (A2C) is logarithmic 



^ = 4 / A ~ |sin(2a)l m^. 
V v ±- b 



Large po/b correspond to x$ ^> 1 where 



(A21) 



1 

Xq = 7T [ U + - 



V± 



2 J y 6A 00 |sin(2a) 



In this limit, 



e n = A oc \sin(2a)\-^e- X0 . 



(A22) 



Thus, for large n and small enough b <§; 7r 2 n 2 £/| sin(2a)|, 
the energy is exponentially close to the gap A^l sin(2a)|. 
The derivative with respect to the impact parameter is 
negative 

de n v±x _ X(j 
db ~ 4b 2 6 ■ 
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